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Abstract 

We characterize the mechanical quahty factor of micro-oscillators covered by a highly reflective 
coating. We test an approach to the reduction of mechanical losses, that consists in limiting the 
size of the coated area to reduce the strain and the consequent energy loss in this highly dissipative 
component. Moreover, a mechanical isolation stage is incorporated in the device. The results are 
discussed on the basis of an analysis of homogeneous and non-homogeneous losses in the device 
and validated by a set of Finite-Element models. The contributions of thermoelastic dissipation 
and coating losses are separated and the measured quality factors are found in agreement with 
the calculated values, while the absence of unmodeled losses confirms that the isolation element 
integrated in the device efficiently uncouples the dynamics of the mirror from the support system. 
Also the resonant frequencies evaluated by Finite-Element models are in good agreement with the 
experimental data, and allow the estimation of the Young modulus of the coating. The models 
that we have developed and validated are important for the design of oscillating micro-mirrors with 
high quality factor and, consequently, low thermal noise. Such devices are useful in general for 
high sensitivity sensors, and in particular for experiments of quantum opto-mechanics. 



INTRODUCTION 

Micro opto-mechanics is earning an increasing interest both for its wide-range applica- 
tions, including high sensitivity measurements of position, acceleration, force, mass, and 
for fundamental research. Among the several breakthroughs obtained in the last few years 



we mention, e.g., radiation-pressure cooling of micro- mochanical oscillators [1|- 



Ml, 



strong 



coupling between mechanical and electromagnetic field variables |7|, \^, opto- mechanically 
induced transparency |9|, [lO|, and, more recently, the observation of the quantum motion 



12|, 



of a nano-mechanical oscillator optically cooled down to its quantum ground state hjj. 

nn 

closely following similar results previously obtained in the microwave region [13|, [1^. Excel 



-ll8| 



lent reviews of this field are published, e.g., in Refs. [15- 

At present, however, quantum effects have been just observed on the variables of the 
mechanical oscillator, while important, long-seek quantum properties of light and of the 
measurement process itself still elude observation. Among them, ponderomotive squeezing 



19-22], quantum correlations, entanglement and quantum non-demolition measurements 
23l-l26|. back-action evasion and measurements reaching and surpassing the standard quan- 



tum limit 



22 
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29|. On such crucial phenomena, only classical simulations are indeed 
reported in the literature [30, |3l| , besides a recent observation of ponderomotive squeezing 
induced by a cloud of cold atoms [32|]. The main reason must be found in the overwhelm- 
ing effects of classical noise sources of thermal origin with respect to the weak quantum 
fluctuations of the radiation pressure. 

Most of the conceived experimental schemes cannot take advantage from laser cooling 
of the main mechanical mode involved in the measurement, while the temperature of the 
thermal bath is the crucial parameter. For instance, the spectral density of the thermal noise 
force is proportional to ks T ^"q^ , where ks is the Boltzmann's constant, T the thermal 
bath temperature, rrif. the oscillator effective mass, aJo/2vr its resonant frequency and Q its 
quality factor. Therefore, experiments aiming to observe quantum properties of the opto- 
mechanical interaction benefit from low mass and frequency (on the contrary, for studying 
the ground state of the oscillator high frequencies are favored by the requirement Huq > 
ksT). Moreover, we remark that the mechanical quality factor plays a crucial role, and 
that experiments take advantage from high optical quality that allows to build high Finesse 
cavities. Finally, high levels of input laser power give stronger radiation pressure, whose 



fluctuations more easily prevail over thermal noise, provided that the oscillator can be kept 
at low temperature. 

Several kinds of oscillators are currently experimented, such as breathing whispering 



galleries 
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9|], thin membranes working as refractive oscillators within high-Finesse cavities 
34| . photonic crystal opto-mechanical cavities |lll. Il2|. Focusing our interest to Fabry- 



Perot interferometers with oscillating mirrors, we mention that very low mass oscillators have 
recently been conceived and tested, based on free-standirig dielectric multi-layer reflectors 



36|, or on tiny mirrors on shaped thin membranes 
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In this article we explore a different approach, focusing on thicker silicon oscillators with 
high reflectivity coating 
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39| . The higher mass is compensated by the possibility to 



manage high power at low temperatures, thanks to the favorable geometric factor 



thicker 



connectors) and the high thermal conductivity of silicon at cryogenic temperature 40|]. With 



such design, it is important to control the mechanical dissipation in the coating layers, in or- 
der to exploit the potentially very high Q of silicon oscillators 
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42| . Here we report on the 



losses measured in a micro oscillator covered by a highly reflective coating, and discriminate, 
by the use of Finite Element analysis (FEM), the contribution of thermoelastic dissipation 
and coating losses. In particular, we test an approach to the reduction of mechanical losses, 
where the size of the coated area is reduced as much as possible to reduce the strain in this 
highly dissipative component and the consequent energy dissipation. 
The loss angle of an oscillator is deflned as: 

, AWt 



2nW 



(1) 



where AWt denotes the energy dissipated per cycle of vibration and W denotes the strain 
energy stored in the elastic body. Mechanical losses are due to the coupling between the 
normal modes and the thermal bath. Several types of losses can be identifled, each one as- 
sociated with a speciflc dissipation mechanism. The measured loss is the sum of dissipation 
arising from different sources Sj, so that we can write AWt = Xlj ^^^i- Note that each kind 
of coupling contributes to the random excitation of the oscillator by drawing from the disor- 
dered motion of the bath, called thermal noise. Until recently, experimenters have focused 
their attention almost exclusively on homogeneous losses, i.e., losses that are described by 
an imaginary part of the Young's modulus in a homogeneous body. This is the case for 
dissipations induced by homogeneously distributed impurities and dislocations. Tools for 



the evaluation of inhomogeneous distributed losses and thermoelastic losses are now avail- 



able, based on analytical models or Finite Elements techniques 43|-|45| . These developments 



allow an accurate estimation of the quality factor of real experimental devices. 

SETUP AND CHARACTERIZATION 

Micro-resonator design and fabrication 

Our micro oscillators are fabricated in a silicon-on-insulator (SOI) wafer where the handle 
layer is 500 um while the device layer is 70 um with a buried oxide of 2 /im. An example 
is shown in Fig. [H Each oscillator is composed of a central round mass suspended by 
eight beams. The oscillator is obtained by engraving the handle silicon layer on one side 
and patterning the structure in the other side by a Deep-Reacting-Ion Etching (DRIE) 
process. The buried oxide is used as etch stop layer. An outer wheel, fabricated on the full 
thickness of the wafer, reduces the mechanical coupling between the main oscillator and the 



high-frequency modes of the whole sample 30|, |46|] . Groups of seven double wheel micro 



opto-mechanical oscillators are arranged in 35 x 35 mm sectors on the wafer, with slightly 
different dimensions to cover a frequency range from 160 to 320 kHz, as reported in table 
m The beams connecting the central mass with the wheels are either 50 or 100 yum wide, 
with length ranging from 300 to 450 fim, and the diameter of the oscillating mirror is either 
700 or 1000 ^m. On the front side of the wafer, a deposition of 17 alternate Ta205/Si02 
quarter-wave layers for a total thickness of 5.244 um provides the highly reflective coating 
(deposited at Advanced Thin Films). We consider two different configurations of the coating: 
in the first case the coating covers the whole front face of the wafer, including mirror and 
supporting beams, while in the second case we used a metallic hard mask to coat only the 
central mirror. The hard mask alignment error was about 100 /xm. 

Opto-mechanical sensing 

The resonant frequency of each device is determined by the noise displacement spectrum 
of the mirror's surface. In order to measure it, we realized a Michelson interferometer 
with a balanced homodyne detection scheme (see Fig. [2]). The laser radiation is provided 
by a cw Nd:YAG laser operating at A=1064 nm. After an optical isolator, a polarizing 



beam-splitter (PBSl) divides the beam into two parts, orthogonally polarized, forming the 
Michelson interferometer arms. On the first one (reference beam) an electro-optic modulator 
(EOM) and a piezoelectric-driven mirror Mp are used for phase-locking the interferometer 
in the condition of maximum displacement sensitivity. A double pass through a quarter- 
wave plate rotates by 90" the polarization of the reference beam, which is then refiected by 
PBSl. The polarization of the beam sent to the micro-mirror (sensing beam), is instead 
rotated by a double pass through a Faraday rotator. The sensing beam is focused with a 
waist of 80 /im on the coated oscillator, and after refiection it is thus totally transmitted 
by PBSl, after which it overlaps with the reference beam refiected by Mp. The overlapped 
beams are then monitored by homodyne detection. The homodyne detector consists of a 
half-wave plate, rotating the polarizations by 45°, and a polarizing beam-splitter (PBS2) 
that divides the overlapped beams into two equal parts sent to the photodiodes PDl and 
PD2, whose outputs are then subtracted. The signal obtained is a null-average, sinusoidal 
function of the path difference in the interferometer. Such a scheme is very weakly sensitive 
to laser power fiuctuations. The difference signal is used as error signal in the interferometer 
locking servo-loop. Low-frequency stabilization (up to a few hundred Hz) is carried out by 
the piezo-mirror Mp, while the wide frequency response of the EOM allows us to extend the 
locking bandwidth up to about 30 kHz. The same signal is sent to a digital oscilloscope 
with integrated Fast Fourier Transform analysis, which provides the noise spectrum. The 
measured spectrum Sy, in V^/Hz, is calibrated through the expression Sxx = Sv{\/'27iVpk)'^, 
where Vpk the peak-to-peak value of the interference fringes and Sxx is the displacement noise 
spectrum in m^/Hz. 

The maximum record length of our scope (250 kSamples) allows a maximal resolution 
of 10 Hz (when sampling at 2.5 MSamples/s to avoid aliasing). In order to analyze the 
narrower resonance peaks, we have also used a digital lock-in amplifier, whose internal local 
oscillator was tuned at 110 Hz from the interested peaks. The beat note, filtered by the 
output integrator of the lock-in with a time constant of 640/is, was then analyzed by the 
scope with a resolution of 0.1 Hz. 

Near the resonance of interest the resonator transfer function Tx{uj) is well approximated 
by the harmonic response 

^-(^) = — r-2 — ^^—^ T (2) 

rrie [00 — oOq + icpt ujooq\ 
6 



where (pt is the total loss factor. The associated noise spectrum of thermal origin at the 
temperature T, is given by the fluctuation dissipation theorem as: 

and its experimental recording can be used to estimate Uq, rrie and (pf The value of rrie is 
then compared with a Finite-Element-Modeling (FEM) estimate to validate the hypothesis 
of thermal origin for the displacement spectrum. 



Experimental results 

In figure [3] (upper panel) is shown the experimental displacement spectrum for one of 
the devices (Dl with homogeneous coating). The resonance peak is well above the noise 
floor and allow the estimation of the relevant parameters. The flgure also report the fltting 
curve according to Eq. ([3]) (plus a flat background). As a comparison, we also show in the 
flgure (lower panel) the spectrum obtained from a similar oscillator, but lacking of the outer 



isolation wheel. As already observed for membrane oscillators 46|, the coupling with the 
background silicon modes splits the main oscillator mode into several peaks, increasing the 
effective mass and reducing the mechanical quality factor. 

In figure m we display the experimental and simulated spectra for a micro- mirror resonat- 
ing at about 215 kHz (device D4). The two experimental curves refer to devices with full 
coating (i.e., coating over the mirror and the beams; FC) and central coating (i.e., coating 
just over the mirror; CC). The resonant frequency in the full coating case is higher, as the 
coating layer contributes to the stiffness of the supporting beams. This effect can also be 
seen in the spectral curves simulated by our Finite- Element model (described below). The 
effective mass rrie is evaluated from experimental data as (72 ±5) fig and (79 ±5) jUg, respec- 
tively in case of the central coating and full coating. These value are in good agreement with 
the FEM estimates (79 ± 5) fig and (83 ± 5) fig, confirming the thermal origin of the noise 
spectra. As a comparison we note that the physical mass of this mirror can be estimated 
by its size as 128 fig; the effective mass is lower because the mirror is not moving as a rigid 
body, but it bends during the oscillatory motion as shown for instance in figure El In table 
H] we report the parameters measured from experimental data for ten devices. We note that 
the quality factor of each CC oscillators is significantly higher than that of the corresponding 



FC version. We shall discuss in Section „ the reason of this behavior and the origin of the 
observed losses in our devices. 

MECHANICAL LOSSES 

In the model of Eq. |2] the measured loss factor accounts for all of the dissipative phe- 
nomena active in the resonator, as defined in Eq. [T] In our case the relevant phenomena 
are thermoelastic loss and structural loss of the various components of the device, as energy 
losses through the support are negligible thanks to the isolation wheel. We consider in the 
following only losses from silicon and optical coating, as the buried silicon dioxide layer is 
removed during the process from the backside of the mirror and of the beams. 

Thermoelastic loss 

Thermoelastic dissipation was first investigated by Zener J47| : in the presence of a nonzero 
coefficient of thermal expansion, when a solid undergoes a vibration other than pure tor- 
sion, the strain field generates a thermal gradient and thus a heat fiow which dissipates 
elastic energy. This fundamental mechanism sets the loss in MOMS devices and precision 
instrumentation at room temperature, for this reason it is the subject of an active area of 



experimental 48|, |49|, theoretical 50J and numerical J45| research. Only in the case of pure 
fiexure the loss factor can be calculated analytically as (j)z = '^c '^ i+tlv^ ' w^^re a is the 
thermal expansion coefficient, Cy the specific heat per unit volume of the material, Y the 
Young modulus, p the density and T the temperature. The oscillator thickness h is involved 
through the thermal relaxation time tz = — |-^, where k is the thermal conductivity. Even 
if the real losses depend on the geometry and the anisotropy of the elastic structure and must 
be evaluated by FEM, these equations give some insight on the behavior of the resonator. 
For instance, in the case of a silicon cantilever with h = 300/im, we have tz — 5.6 /xs and the 
expected loss angle at 250 kHz is about 2 x 10~^. This figure limits the Q factor of a silicon 
fiexure to Q < 4 x 10^ at room temperature, while better performances could be achieved at 



cryogenic temperatures thanks to the changes in the thermal properties of the material 51 |. 



We also note that at these frequencies the thermoelastic dissipation of the optical coating 



is negligible, as its average thermal conductivity is 10 times smaller than in silicon 52| and 

8 



the thermoelastic heat flow is accordingly smaller. 



Structural loss 



The structural dissipation is directly related to the imaginary part of the Young's modulus 
in a homogeneous body, but it cannot be easily evaluated if the device is made of parts with 
different structural losses. If 0s (r) is the loss factor at the position r, the energy dW 
dissipated in one cycle in the volume element dV is simply E{r)(f)s{r)dV, where E{r) is the 
energy stored in the volume element during the motion. In the device the total dissipated 
energy on one cycle is 

AWs= f E{r)<j),{r)dV , (4) 

where both the energy density and the loss factor depend on the position. As a consequence 
the total loss depends on the shape of the displacement within the resonator: modal shapes 
involving large strain in more dissipative parts imply higher losses than modal shapes where 
the same parts are less strained. 

Finite-Element models 

From the above discussion, it is clear that the dissipative contribution from different 
phenomena can be identified only by the comparison with simulations produced by a FEM 
software (in our case, ANSYS Multiphysics) . The Finite Element model of the device is 
based on a three-dimensional 20-node solid element, and the mechanical response of the 
mirror+beams structure is evaluated when it is driven by an harmonic pressure over the 
surface of the mirror. The coupled-field thermoelastic analysis is repeated in the frequency 
range of interest and the convergence of the results is checked against the mesh density. In 
a harmonic thermoelastic analysis, the energy dW lost per cycle in an element with volume 



dV is related to the irreversible entropy produced per cycle at the temperature Tq J45|: 



dV r 36 ^ 39 ^ 
dWte = ^f -JT- kij T^dt (5) 

-'O J cycle OXi CXj 

where 9 is the temperature field, kij the thermal conductivity tensor and Xi the spatial co- 
ordinates; repeated indexes are summed over. This quantity was evaluated for each element 
of the model and integrated over the volume of the resonator to obtain the expected loss 



factor (pte = AWje/(27rl4^). We point out that, to be consistent with the measurements, the 
spatial distribution of the apphed pressure is the same as the laser beam intensity profile 
(Gaussian shape with a waist of 80 /im), and the resulting displacement of the mirror's 
surface is weighted by the same Gaussian profile. The model can simulate, separately or in 
a cumulative way, both kinds of dissipation under study. The thermoelastic loss is evaluated 
from the material properties with no free parameters. On the contrary the evaluation of 
structural dissipation requires as input the loss angles of silicon and of the optical coating. 
The dissipation of silicon at room temperature is strongly dependent on the size of the 
sample, mainly due to thermoelastic losses and surface effects. When the thermoelastic loss 
can be made negligible, the underlying structural loss of silicon wafers is found well below a 



value of 6 Si = 10 
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42| . that we take as reference in our calculations. We point out that 
our results will be almost independent from this number, as losses from other parts of the 
resonator are orders of magnitude higher. Indeed loss angles in the range = (3 — 6) x 10~^ 
53l . |54J | represent the state-of-the-art for optical coatings like the one on top of our silicon 
surface, in spite of the large amount of theoretical and experimental developments carried 
on by the scientific community interested in gravitational wave detectors. 

For each configuration of the coating layer, the loss angle is obtained from FEM data 
using the following procedure: 

• estimation of the resonant frequency of the mode under study; 

• simulation of the transfer function with a given homogeneous structural loss ((/)( = 
10^^) and evaluation of resonant frequency uq and effective mass me from the imagi- 
nary part of the transfer function, according to equation [3l 

• simulation of the transfer function with inhomogeneous losses and evaluation of the 
resulting loss angle 0s : FEM data are fitted with the curve reported in equation |3l 
with given cuq and nie and the loss angle as the free parameter; 

• simulation of the transfer function with thermoelastic loss and evaluation of the re- 
sulting loss angle (pte- FEM data are fitted with the curve reported in equation [3l with 
given uo and me and the loss angle as the free parameter. 

In these simulations we used standard values for the mechanical parameters of silicon and 
the value of Ye ^ 74 GPa for the Young modulus of the coating, measured as described in 
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Section „ We used for the loss angle the value 0c = 6 x 10 ^, within the range of t 



le values 



54| 



measured for this kind of coating in a few experiments at room temperatures [53|, 

In figure [5] we show the simulated noise spectra for a single device (D8 full coating) with 
different kinds of losses. The agreement with the Eq. |3] is generally good and we estimate 
in about ±5% the uncertainty in the evaluation of the loss angle. This error is mainly 
determined by the accuracy limit in the measurements of the dimensions of the device and 
by the discretization errors of the Finite Elements model. The total loss angle is obtained 
as (pt = 4>s + (pte and the expected quality factor is Q = 1/0*. 

To give a better insight on the dissipative behavior of different parts of the device, we 
plot over the displacement the density of the energy dissipated in an oscillation cycle by 
thermoelastic and structural losses. In the case of thermoelastic loss we plot the dissipated 
energy AWte from Eq. [5] (further integrated over the volume), while in the case of structural 
loss we plot AWg given by Eq. |H In figure [6] we display the behavior of the oscillator in the 
case of full coating, while figure [7] highlights the energy dissipated in a micro-mirror with 
central coating. These plots show at a glance that thermoelastic loss is mainly due to the 
flexure of the supporting beams and of the mirror itself, while structural loss is due mainly 
to dissipation in the coating layer. For this reason losses are expected to be significantly 
smaller in the devices with the central coating. 

Comparison w^ith experimental results 

In tableiniwe show the results of the simulations compared with the experimental data for 
a few devices resonating at different frequencies spanning over our full range (165-325) kHz. 
The agreement between simulation and data is within 20% in the case of central coating and 
about 30% in the case of full coating. It is remarkable that such agreement is obtained by 
adjusting a single parameter, namely the loss angle of the coating, within its typical range. 
On the contrary, in the literature, the numerical agreement between simulations and data in 
dissipative micro-systems is usually obtained by fitting the data with a constant amount of 



energy loss, that results to be either of just the same magnitude [55!] or even about 10 times 
bigger J36| than the simulated losses. This extra-loss accounts for some unmodeled energy 
leakage through the support system. 

Therefore we can say that our results confirm the validity of our model based on thermoe- 
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lastic and structural inhomogeneous losses. Our main results can be summarized as follows: 

a) the outer isolation wheel efficiently uncouples the dynamics of the mirror from the sup- 
port system; 

b) the reduction of the coated surface allows to reduce the overall loss in the device; 

c) in the case of central coating the loss is equally contributed by the thermoelastic dissipa- 
tion in the silicon structure and by the coating's structural loss. 

EVALUATION OF THE YOUNG MODULUS OF THE COATING 

The availability of the measurements on two sets of identical devices, only differing in the 
extension of the surface covered by the coating, suggests a way to estimate some structural 
properties of the coating itself. Indeed, the coating layer do not only contributes to determine 
the energy loss as described in section „ but also, through its mass and its elastic modulus, 
it affects the dynamical properties of the resonator. Specifically, in the oscillators with full 
coating the deposited oxide layers cover the beams supporting the mirror, thus increasing 
their stiffness with a negligible increase in the total moving mass, which is mainly determined 
by the mirror. For this reason the FC devices resonate at a frequency higher than CC devices 
(table H]), as we indeed observed in the simulations (an example is shown in figure H]). 

We point out that the oscillation frequency of each device is a function of all of its 
structural and geometric parameters, therefore in principle it could be used to obtain the 
structural parameters of the coating. Unfortunately a typical FEM procedure can evaluate 
the frequencies of the normal modes with a reproducibility not better that 2%, due to 
systematic errors related to the mesh choice and possible biases of the numerical method 
used by the FEM engine. In our case an error of 2% in the resonant frequency corresponds 
to a few kHz, a figure sometimes larger than the observed frequency shift. 

On the contrary the difference Au = Vpc — Vcc is weakly dependent on the choice of 
the mesh and the numerical method, provided that they remain the same in the evaluation 
of both vpc and vcc- Therefore the frequency shift Az/ is mainly dependent on the Young 
modulus of the coating and on possible thickness differences among the two sets of devices 
(FC and CC). We can assume that the actual values of density and Young modulus of the 
coating layer are the same for all of the devices, as they were produced in a single lot. 
Also the dimensions of the devices in the wafer plane are well determined, as they were 
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obtained by a photo-lithographic process followed by a DRIE etching. On the other hand, 
the thickness of each wafer depends on the cutting and polishing processes performed on 
each SOI wafer during the production. 

These considerations can be hardly described analytically, as our devices are not simple 
geometrical structures, but by FEM simulations we can find for each device phenomenolog- 
ical relations between the frequency shift and the relevant parameters. If Yc is the Young 
modulus of the coating. Ah = hpc — hcc with hpc, hcc the thickness of the PC and CC 
device layers, we can write 

Az/ = KyY^ + KtAh + Km (6) 

where Ky-, Kh, Km are constants, measured respectively in Hz/GPa, Hz/yum and Hz, ob- 
tained by the fit of a number of FEM simulations. Here Ky models the sensitivity of Az/ to 
changes of the Young modulus of the coating, K^ to changes in the thickness of the device 
and Km describes the effect of the mass of the coating layer covering the supporting beams. 
The actual values of the constants depend on the geometry of the specific device, therefore 
we have a set of ten equations like Eq. ([6]) (one for each device considered). If the thickness 
of the silicon device layers are known, from these equations we can estimate the Young 
modulus of the coating. 

The K constants that we have found with FEM are reported in table IIIII Our Finite- 
Element models use for the coating a density pc = 4525 kg/m^ and a Poisson modulus 



(Jr. =0.2. The density pc has been calculated from the densities of the Ta205 |56| and Si02 



57l | layers, given that our coating is made of 17 pairs of Ta205/Si02 layers, which amount 
to 2184 nm of Ta205 and 3060 nm of Si02. 

From our data it is clear that, as Kh is on average ~ 2 kHz//im and the frequency shift 
is on average ~ 4 kHz, the thickness difference Ah must be measured with an uncertainty 
better than 0.5 pm to allow a 25% estimate of Az/. Therefore the thickness of the device 
layers hpc and hcc must be measured with a precision better than 1%. The thickness was 
measured on samples obtained by of the actual wafers used for the production of each series 
of devices. Two samples (one PC and one CC) having dimensions 1 cm x 1 cm were cut 
with a quartz saw from the SOI wafers and then clamped and aligned in the same holder in 
order to avoid measurement errors due to the reproducibility of the positioning in the SEM's 
holder. We have used a Jeol JSM-7401P SEM with a magnification of 1500x to have a full 
view of the device layer and the buried oxide. In figure |8] the two samples are shown. A 
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white region where electrons are scattered due to the oxide layer is shown on the left side of 
the two images. All pictures were taken at a the same tilt angle of the electron beam of 3 ° 
(to reach the same experimental conditions). The thickness of the samples were evaluated 
as an average over measurements performed from the top to the bottom of the images. In 
this way we could reduce the effect of the surface irregularity near the interface between the 
device layer and the oxide, which do not give thickness variations greater than 0.2 /xm. We 
have obtained hpc = 70.1 ± 0.2/im and hcc = 69.75 ± O.lyum. 

Once the evaluated K constants and the measured A/i are inserted in equation 0, we 
obtain a set of ten equations for the shift Av as a function of Y . A comparison with 
the experimental data, using a best fit approach, allows us to finally obtain for the Young 
modulus of the coating l^c = 74 ± 20 GPa. In table IIIII we also report for each device 
the estimated Az/jg^, obtained with the best fit of Yf., compared with the corresponding 
measured values Az/gxp- 

CONCLUSIONS 

We have measured the mechanical losses in a micro oscillator covered by a highly refiective 
coating made of 17 pairs of Ta205/Si02 layers. These devices have been developed to detect 
radiation pressure coupling between a low-mass moving mirror and an incident light field 
in a Fabry-Perot cavity, an application where it is crucial to reduce as much as possible 
the thermal noise. Our approach exploits silicon micro-oscillators built in the device layer 
of a SOI wafer. Their design include an embedded isolation stage to limit the mechanical 
coupling to the background wafer. With respect to other kinds of oscillating mirrors, recently 
proposed and tested, that use much thinner structures for the mirror suspension and even 
the oscillating mass, our devices suffer from a higher effective mass. This drawback is 
compensated by the possibility to use high laser power at low temperature, thanks to both 
the favorable geometric factor (thicker links) and the high thermal conductivity of silicon at 



cryogenic temperature [40|. For instance, at 4.5 K our oscillators can bear several mW of 
absorbed power keeping the temperature spread within the device below 1 K. Considering 
that typical absorption in high reflectivity coatings is of few ppm, the intracavity laser power 
can reach ~ 1 kW (obtainable with the typical Finesse of ~ 30000 and an input power of 
~ 100 niW). We remark that temperature homogeneity is important to avoid effects of 
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non-equilibrium thermal noise 



The devices that we have presented are also characterized by a good sturdiness, useful 
for possible applications to high sensitivity sensors. For instance, we can clear them with 
standard cotton-tipped sticks to routinely obtain optical cavities with Finesse around 40000, 
using a standard 100 ppm transmission input mirror. 

In order to optimize the efficiency of the micro-mechanical oscillators in quantum optics 
experiments, it is important to design devices with reduced mechanical dissipation. In this 
work we analyze in particular the advantage brought by a reduced coating area. The results 
are discussed on the basis of an accurate analysis of the possible losses and validated by 
a number of FEM simulations. Thanks to the use of two families of devices (full coating 
and central coating), the contribution of thermoelastic dissipation and coating losses could 
be separated and the observed change in quality factor is in agreement with the expected 
value, demonstrating that the reduction of the coated surface allows to reduce the overall 
loss in the device. We also note that the absence of unmodeled losses confirms that the outer 
isolation wheel efficiently uncouples the dynamics of the mirror from the support system. 

The resonant frequencies evaluated by FEM are also in good agreement with the experi- 
mental data, within the reproducibility observed by measuring different samples of the same 
device. By comparing the resonant frequency of full coating and central coating devices, 
we could estimate of the Young modulus of the coating, in agreement with measures per- 
formed on completely different systems. Indeed, the coating applied over our mirrors is the 
same currently used in optical cavities for metrology experiments and in the large mirrors 
of gravitational wave detectors. As a final remark we observe that thermoelastic loss could 
be reduced well below 10^^ by cooling the system at liquid helium temperatures |51|, where 
therefore our devices should reach a quality factor exceeding 3x10^, at the state of the art 
for micro-mirrors in our frequency range. On the other hand our design could be further 
improved, as most of the losses at low temperature occur in the coating covering the mirror, 
due to the bending of the mirror itself (see figure [7]). A design where only the beams bend 
during the oscillation, while the mirror moves as a solid body with negligible strain, could 
be very useful to further reduce the energy loss in the coating. The FEM models described 
in this work and validated by the comparison with the experiment represent an essential 



tool for such designing activity. The requirements proposed, e.g., in Ref. [30| for obtaining 
pondero-motive squeezing (in particular, a Q of ~ 10^ in an oscillator with a mass of 50 /xg 
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and a frequency of 100 kHz operating at liquid helium temperature) are achievable by such 
design. 
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TABLE I. Parameters of the different micro-mirrors as obtained by the fit of experimental thermal 
noise spectra with the expression ([3]). The experimental errors are estimated from the reproducibil- 
ity of the data obtained with different samples of the same device. 



Oscillator size 



Full coating 



Central coating 



Central disk Beam Beam Freq. (kHz) me(/ig) Q Freq. (kHz) rrie (/ig) Q 

diameter (^m) width (^m) length (/zm) ± 1% ± 10% ± 10% ± 1% ± 10% ± 10% 



Dl 


1000 


100 


400 


167.7 


82 


6300 


164.7 


76 


16300 


D2 


700 


100 


300 


324.7 


49 


5200 


317.3 


50 


13400 


D3 


700 


100 


350 


279.0 


59 


6000 


273.1 


47 


16200 


D4 


1000 


100 


300 


215.2 


79 


6500 


211.7 


72 


17300 


D5 


700 


100 


400 


242.3 


59 


6200 


238.0 


46 


19700 


D6 


700 


50 


400 


208.1 


45 


5200 


203.3 


47 


17200 


D7 


700 


50 


450 


181.5 


50 


6100 


177.4 


47 


16700 


D8 


1000 


100 


250 


249.4 




5900 


242.1 




15100 


D9 


1000 


100 


200 


288.9 




4100 


280.5 




11900 


DIO 


1000 


100 


150 


338.2 




4900 


328.6 




8200 



TABLE IL Comparison between experimental and simulated Q values. 



Full coating 



Central coating 



^s 0te QfC Qexp 4>s 4>te QcC Qea;p 

± 5% ± 5% ± 10% ± 10% ± 5% ± 5% ± 10% ± 10% 
Dl 6.2x10-^^ 4.5x10-^ 9300 6300 2.3x10"^ 4.5x10"^ 14800 16300 
D2 7.4x10-^^ 2.8x10-^ 9700 5200 3.4x10"^ 2.8x10"^ 16200 13400 
D4 7.1x10-^^ 3.7x10-^ 9200 6500 3.2x10-"^ 3.7x10--'^ 14400 17300 
D5 7.1x10-5 3.4x10-5 9500 6200 3.0x10-^ 3.4x10"^ 15600 19700 
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TABLE III. Proportionality constants K describing the dependence of Av from the relevant pa- 
rameters, according to equation [6l Ke is the sensitivity to changes of the Young modulus of the 
coating, Kh to changes in the thickness of the device and Km describes the effect of the mass of 
the coating layer covering the supporting beams. The constants are evaluated by fits on the results 
of FEM simulations. In the last two columns we report the expected values of the frequency shift 
Az^jem, and the corresponding experimental measurements Ai/^xp- In the evaluation of Ai/f^^ we 
have used Ah = 0.4 it 0.3 /um, directly measured, and the best fit evaluation y^ = 74 it 20 GPa. 



Ke 


Kh 


Km 


Az^/em 


Al^exp 


Hz/GPa 


Hz//im 


Hz 


kHz 


kHz 



Dl 78.6 ± 0.1 2258±3 -2845 ±3 3.8 ±2 3.0 ± 3 

D2 105.8 ± 0.3 4405±6 -1826 ±4 7.6 ±2 7.5 ± 4 

D3 110.2 ± 0.2 3868±6 -2130 ±4 7.5 ±3 6.7 ± 2 

D4 83.5 ± 0.1 3007±4 -2885 ±9 4.4 ±3 3.5 ± 3 

D5 92.2 ± 0.2 3399±5 -2108 ±3 6.0 ±3 6.2 ± 2 

D6 97.2 ± 0.2 2921±3 -1004 ±2 7.3 ±3 6.0 ± 2 

D7 80.6 ± 0.1 2581±3 -1003 ±2 6.0 ±2 4.1 ± 3 

D8 95.7 ± 0.2 3129±5 -3155 ±4 5.1 ±3 7.3 ± 4 

D9 102.0 ± 0.2 3546±6 -3171 ±4 5.7 ±3 8.5 ± 4 

DIO 108.7 ± 0.4 4041±8 -3102 ±5 6.5 ±3 9.5 ± 4 
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FIG. 1. Design of the silicon micro-mirror. The diameter of the central oscillating disk is either 
700 or 1000 fiui, while the outer diameter of the isolation wheel is 3000 /im. The beams connecting 
the central mass with the wheels are either 50 or 100 //m wide, with length ranging from 300 to 450 
/im, depending by the target frequency of the device. The thickness of the central disk is about 
70 fira, while the thickness of the wheel is 500 fim. In the design shown here the optical coating 
layer (in violet) covers only the central disk; all devices have been produced also in a fully coated 
version. 
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MOMS 




Scope 



FIG. 2. Scheme of the experimental apparatus. 01: optical isolator; FR: Faraday rotator; EOM: 
electro-optic modulator; H: half-wave plate; Q: quarter-wave plate; PD: photodiode; PBS: polariz- 
ing beam-splitter; MOMS: micro opto-mechanical systems. Lenses and alignment mirrors are not 
shown in the scheme. 
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FIG. 3. Upper panel: experimental displacement noise spectrum for a micro-mirror homoge- 
neously coated (device Dl). The thick solid line (red) is the fitting curve given by Eq. ^, plus a 
flat background. Lower panel: displacement noise spectrum obtained from a similar oscillator, but 
lacking of the outer isolation wheel. 
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FIG. 4. Experimental spectra for a micro-mirror resonating at about 215 kHz (device D4). The 
curves refer to devices vi^ith full coating (FC) and central coating (CC). The simulated curves are 
in agreement within 10% with the experimental results. The frequency shift Ai^exp and Az^jem is 
due to the different contribution of the coating layer to the structure's stiffness, as discussed in 
Section _ 
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FIG. 5. Simulated noise spectra for a single device (D8 full coating) with different kinds of losses. 
Fit with theoretical spectra from Eq. [3] is shown only for the calibration curve with fixed loss 
(d> = 10~^) and for the curve with both thermoelastic and structural losses. 
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silicon 




FIG. 6. Contour plots of the energy dissipated in an oscillation cycle, plotted over the modal 
shape for the device D8 with full coating: a) thermoelastic loss, b) structural loss. In both cases 
the insets show separately the dissipation in the silicon wafer and in the optical coating layer. The 
magenta spot in the center of the mirror gives the size of the optical waist. 



silicon 




FIG. 7. Contour plots of the energy dissipated in an oscillation cycle, plotted over the modal shape 
for the device D8 with central coating: a) thermoelastic loss, b) structural loss. In both cases the 
insets show separately the dissipation in the silicon wafer and in the optical coating layer. The 
magenta spot in the center of the mirror gives the size of the optical waist. 



25 




FIG. 8. Device layer images, including the buried oxide (white region), for the central coating (left) 
and full coating (right) SOI samples. These images are used for the evaluation of the thickness 
variations. 
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